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Flutter  at  Nonlinear  Ac-rody..  :.le  Forces 
by 

P,  S,  Lands  and  S.  P.  Strelkov 

Vibrations  of  an  aircraft  wing  during  fli  ht  should  be  considered  as  vibrations 
in  a  complex  heterogeneous  vibrational  system  with  distributed  parameters.  But  quite 
often  to  reveal  the  nature  of  the  flutter  phenomenon,  or  outoexcitation  of  elastic 
vibrations  of  a  wing  in  flight,  are  used  simplified  systems |l, 23*  One  of  such  simple 
arrdhgenc-nts  is  replacement  of  the  wing  by  a  thin  rigid  plate, elastically  fastened  so, 
that  it  can  execute  combined  "bending-torsional*  vibrations, consisting  of  progressive 
vertical  vibrations  of  the  plate  as  a  single  integral  and  torsional  vibrations  of  sane 
around  a  certain  axis.  At  a  proper  selection  of  parameters  such  an  arrangement  does 
quito  veil  represent  vibrations  of  the  right  wing  tip  of  no  small  aspect  ratio.  Ana¬ 
lysis  of  phenomena  on  such  a  simplified  system  allows  to  examine  the  basic  laws  gover¬ 
ning  flutter  and  to  obtain  a  dependence  of  the  critical  velocities  upon  the  main  dyna¬ 
mic  parameters  of  the  uinginatural  frequencies  of  tha  wing  and  the  magnitude  of  "b-  ad 
betwen  bending  and  buckling". 

In  a  majority  of  existing  experiments  the  flutter  phenomenon  is  investigated 
in  linear  approximation., which  is  valid  at  small  vibrations  and  at  not  high  flight 
speeds. If  the  3pecd  of  flighjfc  is  many  tines  greater  than  the  speed  of  sound, then  in 
conformity  with  the  theory ^.,3^ the  aerodynamic  "forces  depend  nonlincarly  upon  the  angle 
of  attack  (Ifewton  law)  and  this  nonlinear  dependence  has  to  be  taken  into  considera¬ 
tion. 

At  very  high  sypersonic  speeds  the  tine  ~  *  b/v,  where  b  -  wing  chord,  v-flight 
speed,  will  be  insignif icantly  snail  in  comparison  with  the  period  of  oscillations  T. 
The  ratio  'v/T  is  proportional  to  the  StruWiali  number  and  in  the  case  under  considera¬ 
tion  very  close  to  zero,  and  consequently  it  can  be  a  mourned,  that  it  is  permissible 
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to  on: lay  the  results  of  stationary  theory, obtained  for  a  constant  in  tine  ancle  of 
at  tack,  for  the  cane  of  variable  angle  at  vibrations.  The  basis  of  this  assumption  will 
be  explained  below. 

Analysis  of  wing  uutooxcitat  ion  conditions  with  the  aid  of  the  above  described 
simplified  schene  offers  the  possibility  of  explaining  all  basic  features  of  the  pheno¬ 
menon,  c ounce te  d  with  the  nonlinear  nature  of  lift/ancle  of  attach  dependence. 

Analytical  solution  of  the  mentioned  problem  in  view  of  the  nonlinear  nature  of 
equations  is  difficult  and  possible  only  in  approximation.  , 

In  this  report  are  given  results  of  studying  bending-torsional,  flutter  of  a  wing  with 
•IJev’ton*  lift  with  the  aid  of  an  electronic  modeling  machine  type  MI-7. 

The  obtained  dependence  of  critical  flight  speed,  at  which  flutter  appears,  upon 
wing  parameters  end  initial  perturbations^  is  described..  In  spite  of  the  fact  that 
all  measurements  pertain  only  to  qualitative  description  of  phenomena,  nonetheless  it 
can  be  assumed  that  the  possible  value  of  critical  rate  of  flutter  can  be  determined 
theoretically  approximately  with  the  sane  accuracy,  as  it  is  done  for  small  veloci¬ 
ties  in  the  case  o?  applying  the  Bubnov-Gelerkin  method  ]V|. 

1,  Simplified  arrangement  for  studying  bending- torsional. vibrations  of  cantilever 
secured  wing  is  shown  in  fig.l^guations  of  vibrations  of  such  u  system  have  the  farm 

mK  +  ik  —  msbP-F.  lb  +  cd  —  mzbii  =*  M  (1.1) 

where  h  -  displacement  of  axis  of  rotation  (rigidity  axis),  $-  angle  of  rotation  of 
wing  around  the  rigidity  axis,  n-mss  of  wing,  k-elasticity  coefficient  of  wing  at 
vertical  displacements  (bonding  rigidity),  c-elasticity  coefficient  at  rotation  around 
axis  (torsional  rigidity  of  wing),  I-nonent  of  inertia  relative  to  rigidity  axis, 
b-vir.g  chord, gb  -distance  free*  rigidity  axis  to  the  center  of  gravity,  ?  ind  I-I  aero- 
iy.vxaic  forces  nd  moment  of  these  forces  relative  to  the  rigidity  axis.  In  equations 
(li)  the  coordinates  h  and  <9  arc  road  fro:*,  the  position  of  equilibrium  in  the  absence 
of  uerodyn  u.iic  forces. 


z 
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In  the  I  in/ton  a.  procricntioa  the-  resultant  of  aorodyn  :.;tc  forces 

in  care  of  a  tine  constant  angle  of  attack  v;ao  a": lied  in  th  ;*ea..<.tric  c-.atcr  of  the 
vine  (if  the-  uing  is  pre rented  in  fern  of  a  flat  plate)  and  equals 

•  /  •/./<y.i*-iir*»«is»n^  .  ^.0.,)  (I.ii)  . 

lire  position  of  equilibrium  in  the  absence  of  aerodynamic  fo:*ces9  L-lcnpth  of 
vinc^.q-air  density,  v-speei  of  flight.' 

iiaacnt.of  ^  er-odyna.de  forces:  relative  to  the  rigidity  arris  equals. 

M~sa  i~  ~'ur  4- it'll  fl  '  (|.:|)  • 

•  vher e6jb  -  distance  'freu.- rigidity  arris  to  the  point  of  applying- the  equivalent  ?.• 

At  a  .variable  -»ngle  of  wing  attack  the  nagni tilde  of  the  oerody Sards  force  acting 
at  the  given  point  of  the  ving,  ic  determined  by  instantaneous  local:. angle'  of  attack-.-©'  •, 
vhich  at  plate  vibrations  ’.dll  have  at  the  distance  ,  zb  fr era  leading  edge, the  following 


valuer 


\!  i*6(-  f  a.  -  -  /,]  (i.t, 

' Conse^uent-iy, the  tr til- lift  affecting  the 
wing  during,  vibrations, equals  .. 

•  •  I  ‘  .  •  /  •  •  •  .  * 

F;  L/iors\si.r 0* >iuu  0*</s.,  (1.5) 


' ;  secured  vdng. 


F  \  0*.*sign  0*f/c  '  (!.(>) 


In  ai>  analogous-  manner  can- also.be -presented  the  nbme.nt.-cf  aerc-lyanc-ic  forces  l-I, 
Substituting  the  value  ■&?  frou,(l.4)in  (1,6),  ve. have  . 


F  -  IJb*\ jo  +  4-[«fc(s  +  0,  -4)  :>]|*x  . 

and  analogously 

.  X  Men ’Jo  +  +  °»  -y)  —  •*]}* 

.  (1,7) 

M  -LfcV.^  je  ;  *[M=  -  «•  -  1)  -  *!}*(=  i- «. 

• 

t)* 

;<  Men  jo  J-  ’  £©ft(c a,  —  JLJ  Aj  dz 

(1.8) 

It  is  apparent, tint  at  snail  z 


sign  of  the  expression 
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,o  r ;  [#&;= ! !  )  /;j} 

'determined  '07  the  of  the  vilua 

;• ::  t^K:-  -?)-q  1.1 

digression- (l,3)  changer.’  strji  at 

m 


0M> 


'■  0  ..  1  f, 

r«V * 

*•  ii  .  7  - 


r<it  | 


If 'a  is  beyond  the'  limits  ■  of  iht«;cration,t!)3n  this  .ohan~c  of  sir-ii -needs  not  be 
consic'orfed.  Consequently  tocalculate  integrals  in  (l,7(  and  (1,8)  wo  .shall  .exsinihe 
tr~o  cases.  •  •  .  ’  ’•  , 

1,  Acsuain.q  that  zQ<\  .  l.thsn,  as  •  it .  can  -ba  -easily : seen,  . 

;  F  /./*,-(< >*•:  2SHS)  Htfn  i  •  -.-l-  S9’  • '  )-.%vJ 

m  .  .1/  •  ^  |.  , !_ )o«  ivu/]V;  '  (i.to) 

>••««[» -•.»'•]  •  '  .  . 
ISe.ro <0  *  ■frequency  of-  bendir.r  vibrations  of  the  -wir<£,  '•  S-yaluc /s-ropar  tionsl  to 

dtrukhcli'  nuxber,  1  -diriensiohliccs  coordinate  of  vertical  ’c'ispla’cenent  of  the  win-,  ; 

dashes  designate  differentiation  in  accordance  '.Jith  r'inionsibnloda  tine  y 

r  '  !-•  ;t.  "rf-  :s  -  “?•  >•  .Y.i:  t/c<\.  -■ 

. At' snail-- values  of-  the  i-trubhsli  nuriber,  'uhea  sfe,  1,  th.’ condition  zQ  l.'ix  ful¬ 

filled  cyeryuheres,’  In’ addition,'  in  expressions . (1.10)  and  (1,11)  it  is  possible  to 

discard  all  nenbors,- containing. derivatives,  because  -  they  will  he  i?.'nll:ih  "•  -isca 
•  ;  ...  -  'g  •  ••  .  . 

with  <£hnd'  A  ,  We  then  obtain ’expressions,  ccincidinp  v.i  th  '-'ith  (S)  and  (3). 

2*  At  Z0^l.  In  this',  care  •  the  zone  of  integration  in  (1,7)  and  (1,3)  should  be  di¬ 
vided  .into  two  zone  si  .0  /  a  •  ant1  zQ/.  z  Z^l.Wo  will  then  obtain 

h  -  /Af,c*j[8*  +  iVff.O.V  -  2  a 

+  (i ■  +  v(«*  -2Sfx')!^«u*-- 

-  i  -B#.  <i#  _  J_  _  l )  l«)‘  }-  -V  (r,  -  1  +  •&%)  'JO'  -  ‘iSHx’U 

x  siB«[o  ■{■ so-  «t.  -  4-)  ~  J  (Ml) 
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.  .»/  -  -  J[e*o,  +  +  it,zi)  -  2Sov o*’  -} 

.■•+  zwr(-f  +  *(«.- -!-)  +  («.- •'; 

t(*V  +•  3*— -47  “  "*•*)[*  7.* -•+-*)  + 

4.VOO'  [±  +  s.(*.  -  4*)  +  (*.  ~  4-)*) 

-  2V0r'(=.  —  I  '  AV(o,  -  ~)  .Vx'jj  .  (j  .12) 


4^ 


fnrr&.xoiiitjL. 


-i:pronsions  (1.11)  and  (1.12)  ore  valid  at  relative!;'  printer  Jtrukhal1  numbers. 

.  la  this  report  is  discussed  only  the  first  a  _.  pr  on inat i on,  c  or  r  e  o online  to  snail 
Strukfcali  numbers. Solution  of  the  problem  in  this  case  is  such  simpler, and  for  many 
practically' i.  .portant  conditions  such  approximation  if  perfectly  cufficient.it  allows 
to  .explain  the-  basic : essential  traits,  of  the  phenomenon  and  the  influence  of  verious 
•  parameters. 

'  Wo  will  urite  equations  (lj  in  a  different  for:-,sol7in£  sans  v.lth  respect  to  h 
r.nd<$  and  introducing  new  dinenc ionless  values  • 

r  .:,T.isr.  -/  ’[r  //JK 

.  Tnen  we  will  obtain.  a,«*j  . «rj*s#  :  -  »3.,3)Mi.*0sij:u0 

.  0*:#y#.':  a*a3x  #/<ri(3-- 3„)  mu*  9  mcii  V  (1.13) 

'  t  -.  ft  \ 

.  2.  3a sic  results.^uatioris  (1.4)  have  been  nodded  on  the  MI-7  type  machine. 

Durin.~  the  nodelinc  in  equation'  (1.12 )  hove  been  introduced  snail  nenbc-rs.detf.-r- 
:ajainc  oscillation  danpinGJjjTnc  SDcnitude  of  these  snail  rk.-itbers.as  shorn  by.experi- 
nentsihss  practically  no  effect  on -.the  ca'Gnitude  of  the  critical  speed. 

.-  Toe  examined  system  is  nonlinear,  consequently  the  conditions  of  excitation  depend 
ripen-  the  initial  state  of  the  system.  Pisasurossnt  cf  the  critical  flight  speed  value. 


at  flutter  does  take,  place,  uas  carried  out  at.  initial  conditions  x(0)  =  0  and 


5(0) 


*  °. 

.  Such  initial  conditions  indicate  the  presence  of  an  initial  deviation  of  the 
system  froa  its  stationary  state  x  and  determinable  by  equation 

4*0’  ~  -  aV/**'**®  •  h  ft”  |  ft..,  s'  --- i/Mir'i  (2.1) 

Tiie  mentioned  initial  perturbation  my  oricinntej  fee:  example , or.  account  of 
tha  sharp  c'u.nce  in  airspeed  as  result  of  sudden  uiad  qust3,  atmospheric  turbulence 
etc,  2otoilod  investications  on  the  behavior  of  the  oysters  at  wrious  initial  per 
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turbations  have  not  been  conducted  in  tills  experiment, and  the  problem  was  solved  only 
at  sufficiently  high  initial  jolt3  for  the  purpose  of  ox.iaini ng  the  effect  of  nonli¬ 
nearity.  For  example,  the  case  x(0)  =  Oj  $(o)  =  0  corresponds  to  the  sudden  rise  in 
flight  speed  freu  zero  to  a  valuej  correspond  in,"  to  a  givon  q,  uldch,  of  course,  is 
practically  absolutely  unrealizable.  It  is  therefore  only  naturally  to  expect,  that 
real  values  of  critical  speed  will  be  higher  than  the  ones  obtained  in  thio  experi¬ 
ment  (at  corresponding  wing  parameters). 

The  case  of  snail  initial  deviation  from  the  equilibria  position  of  the  system 
in  the  presence  of  eerodj-nar.de  forces,  i.e.::(0)<~. x°, n3  was  not  thoroughly 

investigated  on  the  model,  because  in  this  cars  equations  (1.13)  can  be  linearised, by  wri* 

tinS  if  «-«Y  y  jt  —  x°  (2.2) 

'.’here  <3“ and  x°  are  deter: ined  from  equations  (2.1)  wldch  are  easily  3olvablo,by  chang¬ 
ing  sin  &  into  flf 

°°  ;  i  +  V£)V  (2.3) 

Substituting  (2.2)  in  (1.13)  and  considering  y  as  3mall  values, we  will  ob¬ 

tain  for  them  linear  equations 

y'  +  a’y  -f  <i*4’3«|-  —  2<yu*  (1  —  *33,)  b*tf 

+ aU*f  +  a*a3y  =  29(J*3(a  — 3,)  (2-4^ 

A  study  of  equations  (2.4)  for  stability  offers  the  following  equation  for  criti¬ 
cal  value 

'•  as,M«'.V»r'1y  I*  s»)vi«20;  |*  :  -5*  (2.."*) 

where  fij  -  determined  by  expression  (2.3)  at  q  3  q* 

Table 
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s1 
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Accurate  calculation  by  formula  (2*5)  ia  quite  cumbersome, but  ftpprcrcimate  evaluations 
can  be  obtained  easily, since  the  v_lue  depends  v:oalcly  upon  q^leoults  of  calculating 
q m  at  alpha  =  7»Os  0*2 5  for  values  s  0  and^  »  0,25  ore  Hated  in  tableau  this 

table  for  tho  purpose  of  compering  are  given  q^n  values  obtained  on  the  model  and 
showing,  tint  the  model  doe3  quite  well  describe  the  oscillations  in  the  system  uiidcr 
question,  lhe  observed  error  bears  a  systematic  nature  and  is  explained ,nost  likely, 
by  inaccurate  adjustment  of  the  nonlinear  block  of  the  model  issuing  the  func¬ 

tion  sin2j&<  ,  and  also  by  tho  approximation  of  the  calculation. 

At  creator  initial  deviations  of  the  system  fra:  the  stationary  state  of  equation 
(1,13)  ho  linearisation  can  be  made,  because  nonlinearity  here  leads  to  essential 
differences  in  the  behavior  of  the  system  os  compared  with  linear  theory, 
flcsulta  of  investigating  in  this  casej 
obtained  on  model  of  various  wing  para¬ 
meter  s.are  given  in  fig.2-6. 

In  fig,2  is  given  tin  dependence  of  cri¬ 
tical  value  of  relative  pressure  q*  (flight 
speed  is  proportional  to  va lue q^  of 
torsional  rigidity  of  the  wing, character¬ 
ized  by  value  xi^,f or  various  positions  Fig,2, ’dependence  of  critical  value  q  upon 

torsional  rigidity  of  ths  wing  at  =7,^  *’ 
of  wing  CG  relative  to  the  rigidity  axis  0  for  various  CG  positional  <y=  O.l.and 

$  =  0,?_5,  Solid  curve 3  correspond  to  initial 
under  conditions  that  tlw  rigidity  axis  conditions  x(0)=0, $  (0)  =  0. laser  curve  fer 

f>  =  0,25;  dotted  curve  -  x(0)  =  x°,  <y(0)  » 
is  placed  in  the  geometrical  center  of  3* }  dot-dasli  curve  (theoretical)  at  x(0)  * 

x°,^  (0)  =n?.J 

the  wing,  i»««6(j  B  0,  Vhon  the  CG  is  dis¬ 
placed  to  the  roar  from  the  rigidity  axis  (cf  =  0,25), in  the  system  is  observed  a  pe¬ 
culiar  ’resonance*, i.e, critical  spood  decreases  in  this  case,when  the  frequency  of 
torsional  oscillations  is  approximately  doubly  higher  tlan  the  frequency  .of  bending 
oscillations  (xi^2),  This  phenomenon  is  connoctod.  with  the  nonlinear  quadratic 
dependence  of  the  value  of  lift  upon  the  angle  of  attack  and  appears  only  at  conci- 
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derablo  deviations  of  the  initial  position  of  the  system  fro «  the  stationary  state. 
'..Tien  the-  CG  draws  closer  to  the  rigidity  arda  (s.c.jy  =  0.1)  ‘resonance*  disappears 
and  critical  speed  increases.  '..Tien  the  CG  is  s 'lifted  in  opposite  direction  of  the  r^- 
ldity  a:d.s  flutter  does  not  cricinate  at  all. 

i  Fiy.3  corresponds  to  the  position  of  rigidity  a:dLs  by  a  distance^  »  O.25  from  the 

center  of  tho  wing  closer  to  its  forward  odce^Vhen  CG  coincides  with  Geometrical  center 
of  the  wing  (6=fJa  =  0.25)  in  the  system  is  also  observed  ‘resonance* .  '.lien  the  CG  draws 
closer  to  the  rigidity  arris  (tf  -  0.2)  the  critical  speed  value  rises  sharp  iy.  and  at 
“ « ~  0.17  flutter  disappears.  Cn  the.  other  hand  .when  CG  is  shifted  in  o.  pwitc  direc¬ 
tion  (fi—  0.3)  critical  srecd  decreases  ana  ‘resonance*  disappears.  The  dependence  of 
critical  speod upon  the  position  of  CG  in  the  case  under  consideration  at  constant 
:-i2  is  civen  in  fig.4. 


Fip.lj.IteP-naance  of  critical  value  q. 
upon  the  CO  position  <jrat  a  -  7*6  -  0.25. 

:d2  »  3  ^ 

A  change  in  parameter  a  (chance  in  magni¬ 
tude  of  radius  square  of  winy  lncrtia)leada 

?i-.3.Dcpendencc  of  critical  value  $»upon  to  die. lacxcnt  of  •resonance*  frequency,  be— 
the  weight  of  the  wing  at  torsion  at  alpha  oauaey 

*  7,^  *0.25  for  various  eooitions  of  CG*  h t  smaller  a  ‘resonance*  occurs  at  yrcatdfc 
0  »  0?2,  0.25}  0.3 

tor si on. 1  frequency  ratios  to  the  bending 
frequency  (fig.*>).  Reson.jr.ee  becomes  move  pronounced  here.  At  greater  a  values 
■res. nance*  disap. ear  entirely  (o.g.a  a  10). 

In  (fig.6jij  civen  tho  dependenco  q  upon  :dL2  in  case  when  tho  rigidity  a::ia  is 
shifted  by  a  distance  «*  0,25  toward  the  trailin™  edge  of  the  wine. 
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vie, o.Dc penuenco  of  critical  value  q* 
of  '.;i  g  rigidity  ill2  u.-oa  torcion  in  case 
6»“  -  0.25.0C  »  7*f  *  0.1  na*»  0.2 

In  this  case  the  critical  pressure  value 


rises  almost  szxcthly  ani  •rcsonur.ee*  is  not 

7iC*5»j:tePca?-ence^of  critical  value  q# 
of  vine  rigidity  ■l3.<  upon  torsion  atfl^  observed. 

»  0.25  '.ndo  »  0.25  for  various  values 

of  parameter  a.  In  this  e:peri:ient  was  also  i;wa?ti£pted 

the  effect  of  fucelace  on  the  mcnitude  of  critical  speed,at  which  flutter  does  take 
place.  In  this  case  was  used  a  wine  uodel.shoun  in  fi£.7. 

3quaticns  of  oscillations  of  such  a  model  have  the  fom  of 

riiii  /.  [h  -  II)  -  nizlill  h 


•..-here  I\j-nass  of  fusealnj,  ^displacement 
of  fuselage  relative  to  the  equilibrium 
position,  -  rigidity  of  suspension. 

The  investigation  was  made  at  I-^/n  =  10, 
3  0.7  *  0  .25.  Ike  r>..:ulta  showed 

thst  in  this  case  the  consideration  of  the 
fuselage  has  a  very  slight  effect  on  the 
mn.-nitude  of  critical  flutter  rate. 


1 0  -  mzhit  M 

Mjl  I.  [II  -  h)  r  I  JI  «• 


(2«) 


?ifj,7 .Simplified  system  of  wiiy;  with  fuse* 

lace. 


She  critical  3  joed  values  with  consider  at. 'on  of  the  fu-' elace  -re  shcftin  in  .'tJaph 
in  fic.p.hy  th  sif.il attf*  0,2|  #  »  7.  other  parens  ter  r  lues  the  chsrvre  in 
critical  s.-eod  with  consider  .tion  of  the  fuse  la  c*  was  fount!  to  ho  even  more  insigni¬ 
ficant. 
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ARL  (ARB) 


OTHER  AGENCIES 


Nr.  Copies 


1 

1 


MAJOR  AIR  COMMANDS  Nr.  Copies 


AFSC 

SC  FDD  1 

DM  25 

TDBTL  5 

TDBDP  2 

SSD  (SSF)  2 

APGC  (PGF)  1 

ESD  (ESI)  1 

RADC  (RAY)  :  1 

AFWL  (WLF)  1 

AFMTC  (MTU)  3 

ASD  (ASYIM)  2 


DIA  9 
AID  2 
OTS  2 
ABC  2 
PWS  1 
NASA  1 
ARMY  (FSTC)  3 
NAVY  3 
NAFEC  1 
RAND  1 
AFCRL  (CRXLR)  1 


\ 


FTD-TT-63-201/1+2 


11 


